Bound States in NRQCD/NRQED and 
the Renormalization Group^ 

Iain W. Stewart 

S ■ Department of Physics, University of California at San Diego, 

O ■ 9500 Gilman Drive, La JoUa, CA 92093 

(N; 

Oh' 
(D ■ 

^ ': Abstract. 

^*^ . The application of renormalization group techniques to bound states in non- 

relativistic QED and QCD is discussed. For QED bound states like Hydrogen and 

I . positronium, the renormalization group allows large logarithms of the velocity, Inv 

^ ' (or equivalently Ina's), to be predicted in a universal and simple way. The series of 

0^ i (alna)'s are shown to terminate after a few terms. For QCD one can systematically 

t^^ I sum infinite series of the form [as Ina^]'', and answer definitively the question "as at 

—J. . what scale?". 

o: 

S I For Coulombic bound states of two fermions, the relevant scales include the mass 

■^ ■ of the fermions, m, their momentum, p ~ mv, and their energy, E ~ mv^. Since 

p ■ f <C 1 it is useful to calculate properties of these bound states in a double expansion 

' ! in f and a, with terms a/v ~ 1 kept to all orders. However, the expansion may 

D I still involve large logarithms of v or a, which appear through factors of \n{p/m), 

i-^ ' \n(E/p), and \n(E/m). In this talk I discuss how the renormalization group can be 

>• ■ used to systematically predict and sum powers of alnt> [1-6]. 

K> ! At a given order in QED, terms involving Inct typically give the largest contribu- 

H I tions, and the precision of experiments make their prediction quite important [7] . 

■ - - ' In Ref. [5] it was shown for the first time that Ina's in the Lamb shift, hyperfine 

splittings, and annihilation decay widths can be predicted with the renormalization 

group. This is in contrast to the usual method of computing these logarithms by 

evaluating matrix elements at the scale m. The calculations are simple enough that 

we can simultaneously treat Hydrogen, muonium (/i+e"), and positronium {e~^e~). 

For the ti system near threshold, the relevant scales are rrit ~ 175 GeV, 

rritv ~ 30 GeV, and rritv'^ ~ 5 GeV, which are all ^ Aqcd and can be treated 

perturbatively. In QCD there is a strong dependence of the coupling on the scale; 

as{mt) is much different from asimtv"^). The renormalization group allows us to 

handle this complication, or equivalently it allows us to systematically sum terms 
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FIGURE 1. a) Paths in the {^ij,f^s) plane for one-stage and two-stage running. b),c),d) Ex- 
amples of the fxu and fxs dependence of the Fcynnian rules. 

(cts In f)'^. Predictions for the running coefficients of the ti potentials and the ti 
production current will be discussed below [2,4]. 

Effective theories for non-relativistic QED and QCD (NRQED and NRQCD) al- 
low the double expansion in v and a to be performed in a simple way [8]. However, 
application of the renormalization group in these theories is complicated by the 
presence of two low energy scales, p and E, which are coupled by the equations of 
motion, E = p^ /{2m). If one attempts to lower the cutoff on the energy to i? ^ A, 
then one can still excite larger momenta p ^ y/mA. Using dimensional regulariza- 
tion one can deal with this coupling of scales by using a velocity renormalization 
group [1], which has a subtraction velocity u rather than the usual subtraction 
momentum /i. Running in one-stage from z/ = 1 to z/ = w simultaneously lowers 
the subtraction point for momenta, ns = ^^v-, to the scale mv, and the subtraction 
point for energy, /x^ = mz/^, to mv"^ . In Fig. la this one-stage approach is con- 
trasted with the alternative two-stage approach where one ffist runs from fi = m 
to mv and then runs from /i = mv to mv'^. 

In Ref. [6] it was shown that for QED bound states, the most obvious method of 
two-stage running fails to reproduce terms involving (Ina)'^ with k > 2. This occurs 
because in the two-stage method the coupling between the energy and momentum 
is ignored. 

To calculate observables at the hard scale m, the way the effective theory is 
formulated does not matter too much as long as it has a consistent power counting 
in V. Below m, only on-shell degrees of freedom are kept in the effective theory, i.e. 
degrees of freedom which fluctuate near their mass shell. For example, a potential 
gluon exchanged between two quarks has energy ~ mv'^ but momentum ~ mv, and 
is therefore far offshell. Instead of a potential gluon the effective theory has a four 
quark operator 

^P = -^y{P, P') /^5 ^P'^ ^P X-p'^ X-p , (1) 

p,p' 



TABLE 1. QED Ina's which follow from the leading order (LO) and 
ncxt-to-lcading order (NLO) anomalous dimensions in Ref. [5]. /i+e^ predictions 
include 1/m^ dependence, (h.f.s is hyperfine splitting and AF/F is the e+e~ decay 
width correction.) 
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** There is a growing consensus that the value in Ref. [9] is correct [16]. 

where ipp (xp) destroys a quark (antiquark) with momentum p and spin and color 
labels are suppressed. V{p, p') is the potential 
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where the Ui{i')^s are running coefficients. In QCD the color singlet and octet 
coefficients have different values. We have absorbed in the f/j's the dependence on 
the fermion masses; only the momentum dependence is important for the power 
counting. For the QQ and QQ potentials in an arbitrary Lie gauge group, the 
matching coefficients at one-loop to order f^ can be found in Ref. [3]. 

The effective Lagrangian has quarks with {E,p) ~ {mv'^^mv) interacting with 
soft gluons which have {E,p) ~ {mv, mv), and ultrasoft gluons which have {E,p) ~ 
[mv'^,mv'^) (see Refs. [1-4]). To implement the velocity renormalization group we 
renormalize the Lagrangian and compute anomalous dimensions. This procedure 
is fairly simple since it can be done treating the Coulomb potential perturbatively. 
In Figs. lb,c,d the yu's which appear in some typical interactions are shown. Fig. lb 
shows a quark interacting with a single ultrasoft gluon. Physically, the fact that 
the ultrasoft mode involves the coupling g{fj,u) makes sense; due to the multipole 
expansion the scale mv'^ is the only scale it sees. Fig. Ic shows a soft gluon scattering 
off a quark, and Fig. Id shows an insertion of the potential. For these interactions 
the parameter ns ~ i^v appears. 



TABLE 2. LO and LL values of cocfRcicnts of the tt potential in QCD. 
Here Ua — SC/a and the values are in units of the top mass mt- 

Ui'^mt Ul^'^m'i U^'^mj U^/^m^ U^Amj U'^'^mj 



v = \ 

V = V 



-^^M =t:m d [tbd [ttb 27nr 

-0.03 -1.49 0.63 0.53 0.16 3.11 



Below the electron mass the electromagnetic coupling in NRQED does not run 
but coefficients in the potential do. For fermions with mass and charge (mi, — e^ 
and (m2, Ze) we find the anomalous dimensions: 
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2 / mim2 ' 

where U2+S = U2 + UsS'^ and the p^'s are mass dependent numbers [5]. Solving these 
equations gives the results summarized in Table 1.^ Taking the matrix element of 
the leading log (LL) value of U2{i') gives the a^lna Lamb shifts for Hydrogen, 
muonium, and positronium. Furthermore, the LO anomalous dimension is inde- 
pendent of V so there are no higher terms, Q;'^"'"^ln a for k >2. At next-to-leading 
log (NLL) order the most logarithms are generated by the Pc22^c^2(^)^ term which 
gives the a^ In^ a Lamb shifts. Thus, there are no terms a^+^ In a for k > A. 
Hyperfine splittings are generated by the S^ terms in Eq. (3), and the ortho and 
para-positronium widths are generated by imaginary terms which enter through 
the matching condition for t/2+s(l)- 

NRQCD is better at generating logarithms than NRQED since the running of 
as causes all potential coefficients to run. For ti the change in the color singlet 
couplings from z/ = 1 to z/ = f = 0.15 are shown in Table 2.^ The largest changes 
occur in the spin independent couplings UJf (u), U^^\iy), and U2 (i^)- It is these 
couplings which depend on as{mtv'^) since their anomalous dimensions have contri- 
butions from ultrasoft diagrams. In Fig. 2a we plot the two-loop running of UJf (u) 
whose value changes by an order of magnitude between u = 1 and u = v. The full 
theory ti production current gets matched onto a current in the effective theory 
iYt = Ciip^cr'^X-p + • • •• At NLL order the running potentials mix into the running 
of the production current coefficient ci(z/) [1]. Fig. 2b plots the running of ci(z/) at 



'^' Note that the a^ln a Lamb shift requires the LL running of V^^' since this potential mixes 
into C/2, and the a^ In a Lamb shift requires the NLL running of Uk- These will be discussed in a 
future publication. 
^) The renormalization group improved static potential is considered in Ref. [17]. 
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FIGURE 2. For ti production near threshold the running of the color singlet and octet l/|k| 
potentials are shown in a), and the NLL value of the production current is given in b) (solid line) 
[4]. In b) the large, medium, and small dashes are the LO, NLO, and NNLO [18] matching results. 
In a) and b), i^ = 1 is the scale n = nit and the solid vertical line is the Coulombic velocity. 

NLL from Ref . [4] . Summing the logarithms improves the convergence by reducing 
the size of the NLO matching coefficient by a factor of 2. It would be interesting 
to see if a similar improvement in the convergence takes place for the rather large 
NNLO matching correction found in Ref. [18]. 
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